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The spin dynamics of radical-triplet system (RTS) has been calculated by using the Lindblad 
formalism within the theory of open quantum system. The single-radical-triplet system (SRTS) is 
considered here for single-qubit quantum gate operations while double-radical-triplet system (DRTS) 
for two-qubit operations. The environment effects taken into account include the spin-lattice relax¬ 
ation of the triplet exciton and radical spin-i, the inter-system crossing process that induces the 
transition from singlet excited state to the triplet ground state, and the rather slow relaxation pro¬ 
cess from the triplet ground state back down to the singlet ground state. These calculations shown 
that the line shape broadening is strongly related to the exchange interaction between triplet and 
exciton, which can be understood as a spontaneous magnetic field created by the triplet renormalises 
the original spin-i electron spin resonance spectra. This work will provide key information about 
the spin dynamics for building optically-controlled molecular quantum gate out of radical-bearing 
molecules. Moreover, this has generated the further theoretical question on how the mixture of 
fermion and boson behaves. 
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I. INTRODUCTION 

Molecular-base electronics attracts much attention due 
to the advantages over conventional semiconductors [lj, 
e.g., ease to engineer molecules, long spin-lattice relax¬ 
ation time, and relatively less restrictive production con¬ 
ditions. In this sense, molecular-based electronics is a 
very promising way to tackle challenges obstructing the 
implementation of spintronics and quantum information 
processing (QIP). However, the manipulation of individ¬ 
ual electron spin and controlling the interactions between 
electron spins in molecular species is difficult and crucial 
for the realization of spintronics and QIP in the molec¬ 
ular level. Time-Resolved Electron Paramagnetic Reso¬ 
nance (TR-EPR) SH3 is a very useful tool for observing 
spin polarization and spin-spin interactions both for the 
molecules in the ground state and the excited state Q. 

Against the above background, the radical-triplet sys¬ 
tem (RTS) (H which contains spin-bearing radicals and 
optical-active molecular species could be a prototype for 
molecular-based spintronics and QIP. From the literature 
®, RTS could contain either two radical spins or two 
radical spins and one photo-active coupler, i.e., single¬ 
radical-triplet system (SRTS) or double-radical-triplet 
system (DRTS). However, more radicals or photo-active 
species could be synthesized. TR-EPR is applied right 
after molecules are excited by laser shining from the 
ground state to the triplet excited state through inter¬ 
system crossing. The interactions between these spins 
and the spin polarization effect can be detected by TR- 
EPR. Compared with a large volume of experimental 
works on RTS there is very rare theoretical work to ex- 
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plain the observed EPR spectra from a point of view of 
the open quantum system, i.e., considering the linear re¬ 
sponse of an open quantum system which interacts with 
the environment when it is perturbed by an external field, 
rather than a closed quantum system although there is 
some theoretical work based on a rather ambiguous the¬ 
ory @,0. 

Therefore, in this paper we will first present a general¬ 
ized model HU for the quantum open system in the limit of 
Markov process, i.e., the system forgets its history much 
faster than its evolution. Following this general model, 
we will describe the time evolution of quantum states 
within RTS. And then we will produce some numerical 
results of TR-EPR simulations, and compare them with 
the experimental results. Furthermore we will base on 
our theoretical model to make some predictions which 
is not present in experiments, especially the spin corre¬ 
lations between two ^-spins. At the end, we will draw 
some conclusions and have some discussions. 


II. LINEAR RESPONSE THEORY FOR THE 
QUANTUM OPEN SYSTEM 

First let’s start with a closed quantum system. Sup¬ 
posing we have a closed quantum system described by 
Hq and perturb this system by using an external field 
A f(t)B, where A is a small factor, the linear response 
of the closed system due to a time-dependent pertur¬ 
bation is described by the well known Kubo’s formula 
Q and the time-dependent and frequency-dependent 
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response functions for an observable A read 


= ^(r)<[B/(0),i/(r)]) (1) 

/ +oo 

dT<t> AS (T)exp(iuJT). (2) 

- OO 


However, for an open quantum system we may have to 
access to a coarse-grained description in which the co¬ 
herent Hamiltonian supplemented by incoherent terms 
representing the interaction with the environment. The 
most general description for the evolution of a quantum 
state of an open quantum system M can be formulated 

by 


f = ^-[H i p]+Y J [A,pAl-\{p ) AlA,}} = tp 


p{t) = e c pi 0), 


( 3 ) 


where the equation in the first line is so-called Lindblad 
equation. Here A M describes the quantum jumps due to 

the system-environment interactions and C is the Lind¬ 
blad super operator. In order to derive formulae for the 
linear response of an open quantum system we should 
take some advantages of Heisenberg representation and 
Interaction representation. But before doing this, we can 
envisage that it is quite ambiguous to define the Hermi- 
tian conjugate of a super operator if we want to use the 
interaction representation. Therefore, first we define a 
scalar product of two operators by 


(i,H) = Tr[itB], 


( 4 ) 


where f ensures the normal linearity conditions for a 
scalar product. This scalar product looks just like a vec¬ 
tor scalar product on the ’extended vector’ representa¬ 
tions of the operators A and B. With the help of this 
definition, we can now define a Hermitian conjugate of a 

super operator by stating that t should satisfy the fol¬ 
lowing equations, 

(A, Eb) = (. bA , b) = (b, bAy. (5) 


A. Heisenberg representation 

With the definition stated in equation [5] the Hermi¬ 
tian conjugate of a super operator corresponds to the 
ordinary matrix Hermitian conjugate of the four-index 
matrix representation as follows 


In particular if we look at the expectation value of any 
Hermitian operator O with a time-independent Liuvillian 

(6) [t) = Tr[d pit)] 

= Tr[de £t p{0)\ 

= id\e h P m 
= ie bt 6\ P m 
= Tr[O(f)p(0)]. 


Therefore, 


O(t) = e tU 6 

dtOit) = boy). 


( 7 ) 

( 8 ) 
(9) 


We can also generalize these to the case where the Li¬ 
uvillian is time-dependent as follows 


pit) = Texp[ / £(s)ds]p(0) (10) 

Jo 

O(t) = Texp[ f £^(s)ds]0 (11) 

Jo 


B. Interaction representation 

As shown in the following, we can see that starting 
from the interaction representation we are able to derive 
formulae for the linear response of an open quantum sys¬ 
tem systematically and easily. Let’s first decompose the 

Liuvillian to two parts, i.e., C = Co + £i, where C\ is 

small in some sense compared to £o. For the moment we 
suppose both of them are time-independent. Then we 
can define an interaction-picture operator by 


Oiit) = e £ o*d, 


( 12 ) 


and then we can derive the density operator in the inter¬ 
action picture 


Tr[Op(t)] = Tr[O z (t)p/(f)] 

= (6,e^V/W) 

= Tr [Op At)]. 


(13) 


Therefore we require 


Pi it) = e £ot pit) = e £ot e £t p(0). (14) 

The equation of motion for pi is now 


{&)kl,ij = ibj,klY 


( 6 ) 


d t pi [t) = Cuit)pi it) 

biit) = e- £ot bb ot 


(15) 

(16) 
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C. Linear response theory for stationary reference 

state 

We first derive the formulae in the case where po is 

a stationary state, i.e., e Cot po = Po- Supposing C± = 
A f(t)B, where A is the small parameter, B is a general 
super operator, and /(f) is an arbitrary function of time. 
Now we imagine the perturbation is turned on gradually 
from f = — oo, and look for the change in the expectation 
value of an observable A. We can evaluate this in the 
interaction representation 

POO 

{A)(t) - (A) 0 = X dr fit — T)Tr[A I (t)B I (t — r)po\ 

J 0 

/ oo 

f{t-T)<j>AB{r)d,T (17) 

-OO 

</>Ab(t) = 6 , ('r)Tr[A/(r)Sp 0 ] = 9{t)TA[E* 

Notice that if the perturbation B is a Hamiltonian, 
so B = -C[B,-] and B t = ^[B, •] = — B, then 
cj> ab ( t ) = j-0(t)([B i (0),A i {t)}) o , which is in agreement 
with Kubo’s formulae. Now we can do the Fourier trans¬ 
formation xab(w) = dr/>AB(r)exp(fwr). By using 

the equation e^° T = v^e Aoi T(v^) T , where Aoi are 
the eigenvalues of Co and are the right and left 
eigenvectors of Co- After some algebra, finally we have 

XAs(w-Me) = 1 -[{B(C\ -*w-e) _1 i,p 0 ) 

-((>4 -iuj - e)~ 1 AB,p 0 )\ 

= ^Tr [A(C\ -iui- e)~ 1 Bp 0 
- bA(C\ — iuj — e)~ l po) (19) 

D. Linear response theory for time varying 
reference system 

If po is not stationary of Co, we can still have the sim¬ 
ilar equation for the response function as follows 

/*oo 

(A)(t) - (A) o = A / dr fit — r)(f>(t — t) 

J o 

= 0(r)Tr[Aj(r)Bpos(f - r)] 

9{t)Ty[ABi(—t )pos(t)], (20) 

where posit) = Pa is the unperturbed density opera¬ 
tor in the Schrodinger picture to time t. If we further as¬ 
sume that f{t) = e et cos(uit). After some algebra similar 
to that in the case for stationary initial density matrix, 
we can have the following frequency-dependent response 
function 

Xab(u + ie,t) = ~Ar[A{C\ + iuj - e)~ l Bp S o{t)}{2\) 

which is in agreement with the equation in ref.@. This 
is also the basis for our simulation of TR-EPR spectra. 


III. SYSTEM DESCRIPTION AND 
THEORETICAL METHODOLOGY 

In most of the TR-EPR experiments we are interested 
in molecules are exposed to the environments, in gen¬ 
eral a liquid solution. The spin-lattice relaxation time 
is very long (ps) at room temperature for the radical 
spins. And the system-environment interaction is almost 
through spin-orbital couplings and phonon effect. The 
inter-system crossing is really fast for some molecules due 
to the combination of spin-orbital coupling and electrical 
coupling between singlet excited state and triplet state. 
However the above two mechanisms lead to very slow re¬ 
laxation back to the ground state; the life time of the 
phosphorescence is about ms even at room temperature 
because the electrical coupling and spin-orbital coupling 
is relatively small. In summary, Markovian description 
which is appropriate for weak-coupling situations should 
be a good approximation for RTS. So the main task is 
to build up the effective spin Hamiltonian and quantum 
jumping operators, solve the resulting Liuville equation, 
and evaluate the spectra according to equation [5TJ How¬ 
ever, before that we should also develop the formalism 
for the linear response of quantum open systems as TR- 
EPR experiments involve both system-environment in¬ 
teraction and the linear response of open systems. 

A. Effective spin Hamiltonian 

The effective spin Hamiltonian for SRTS reads 

^SRTS = 9tPBB ■ S + g r p B B • s + JS • s 

+DS 2 z +E{S 2 x -S 2 v ), (22) 

where the first is Zeeman energy; g t is the g-factor for 
the triplet and g r for radical, the second term is exchange 
interaction between radical and triplet; J is the exchange 
constant, and zero-field splitting D and E for triplet. B 
is magnetic field. S is the triplet spin, and s is A -spin. 
Similarly for DRTS the effective spin Hamiltonian is 

^drts = 9tPBB ■ S + g r pgB ■ (si +s 2 ) + JS • (si + s 2 ) 
+DS 2 + E(S 2 x -S 2 ), (23) 

where Si and Si refer to two radical spins. 

If we further include the laser field into the Hamilto¬ 
nian and adopt a rotating wave approximation (RWA), 
the Hamiltonian reads 

Bsrts = |£)[^sRTs]( f l + • s](es| 

+\gs)[g r pBB ■ s](srs| 

+V(\es){gs\ + |ss)(es|) (24) 

for single-radical-triplet system, and 

Bdrts = + \ es )[9rPBB ■ (Sr + s 2 )](es| 

+\gs)[g r p B B ■ (si + s 2 )](0s| 

+V(\es)(gs\ + |ffs)(es|) 


( 25 ) 
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for double-radical-triplet system. Here gs, es, and t 
means the ground state, excited singlet state, and triplet 
state, respectively. V is the electrical dipole matrix ele¬ 
ment due to the laser field within RWA, and we assume 
that the exchange interaction is negligible when the opti¬ 
cal active molecule is in either ground or excited singlet 
state 0. The dominant electrical transitional process of 
TR-EPR experiments for SRTS and DRTS is illustrated 
in Figure ©• 



(a) 



FIG. 1: The schematics of energy levels for different spin 
configurations, the photo-excitation and inter-system crossing 
processes. For each figure the ground and excited singlets of 
optical active species are shown on the left hand side while 
triplets on the right hand side. SRTS is shown in (a) where 
Soa is the spin of optical active species, and Sr refer to the 
radical spin. Similar symbols apply in (b) for DRTS. 


B. Quantum jumping operators 

We are considering two kinds of quantum jumps op¬ 
erators; one is the relaxation process and the other is 
inter-system crossing. The relaxation processes include 
the relaxation of radical spins and the triplet back to 
the equilibrium. The electronic transitions from the ex¬ 
cited singlet state (Si) to the lowest triplet state Ti will 
happen due to the inter-system crossing as shown in Fig¬ 
ure ©> and the electronic transition from triplet down to 
the ground state (another inter-system crossing) is also 
shown in Figure © which brings the system back to the 
ground singlet state. 

Combining the effective spin Hamiltonian and quan¬ 
tum jumping operators, now we can use Liuville equation 
to have the time evolution of the whole system under the 
interaction with the environment. As the Hamiltonian 
part of Liuvillian is just its commutator with density op¬ 
erator, here we list the Liuvillians describing the quan¬ 
tum jumps. Those Liuvillian super operators constructed 
by Lindblad operators are following, 

n * 1 

= ( 26 ) 
/j=i 

For SRTS, i = 1 to 4, and 7 ^ describes how fast the 

quantum jumps happen. C\ with m = 3 is a super 
operator describing the relaxation of 5 -spin; l\ = s + , 
If = s~, and l\ = s z , where s\ z operator is responsi¬ 
ble for the pure dephasing of the quantum states. £2 
with «2 = 3 is a super operator describing triplet spin 
relaxation; If = S + , If = S~, and If = S z . Here we as¬ 
sume that the inter-system crossing due to the spin-spin 
coupling is much faster than that due to the spin-orbit 
coupling. So the total spin angular moment should be 
conserved. Therefore we have the following operators for 

£3 with 11,3 = 2 is a super operator describing the elec¬ 
tronic transition from excited singlet to the lowest triplet 

state; if = ||, |) t (|, \\ es and if = ||,-§)*(§, -§|e»- £4 
with 714 = 2 is a super operator describing the electronic 
transition from the lowest triplet to the ground state; 
It = \\,\)ga(\,\\t and if = \\,-\) ga (\,-\\u where 
spin states are written as the form |s,m s ). 

Similarly for DRTS, i = 1 to 5 as we got another set 
of operators for the relaxation process of the addition 

radical. £1 with n 1 = 3 is a super operator describ¬ 
ing the relaxation of 5 -spin si; l\ = sf, If = and 

^3 = siz- £2 with ri 2 = 3 is a super operator describ¬ 
ing the relaxation of 5 -spin S 2 ; If = sj, If = s 2 , and 

I3 = S2z- £3 with ?i3 = 3 is a super operator describing 
triplet spin relaxation; l\ = S + , if = S~, and I 3 = S z . 

£4 with 714 = 7 is a super operator describing the elec¬ 
tronic transition from excited singlet to the lowest triplet 
state; | 1 , 1 )J< 1 , l| es , [ 1 , 0 )J( 1 , 0 | es , | 1 ,- 1 >J< 1 , — 1 |„, 
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|1,1}?(MU , |1, 0 } t 2 (l, 0 | es , |1,-1)2(1,-1|„, and 

| 0 , 0 ) t ( 0 , 0 | es , where we use superscripts 1,2 as there 
are two triplet manifolds when S oa = 1 ; one is 

from 1 © 1 and another from 1 © 0. £5 with 

11,5 = 7 is a super operator describing the elec¬ 

tronic transition from the lowest triplet to the ground 
state; | 1 , l) gs (l, l\j , | 1 , 0 ) gs (l, 0 | t \ | 1 ,-l) fla (l,-l|J, 
|1> l)ffs(l) 1 It , |l, 0 ) flfl (l, 0 | 2 , |l,-l) fls (l,-l| 2 , and 

|0,0) ss (0,0| t . 

IV. NUMERICAL SIMULATIONS FOR EPR 

EXPERIMENTS OF SRTS AND DRTS 

By now we have Hamiltonians including equation (1241) 
and equation (1251) for the coherent process driven by zero- 
field splitting, exchange interaction, magnetic field, and 
laser field. In addition, the dissipative process of the 
whole system is described by the Livillians defined in 
equation (l26l) . So next we are going to put in some rea¬ 
sonable parameters and calculate the time evolution of 
the density operator, especially the populations for the 
excited singlet state and triplet state and the field de¬ 
pendency of EPR spectra due to a perturbation of trans¬ 
verse microwave field, i.e., a linear approximation for 
the EPR spectra. In general, the relaxation of triplet is 
much faster than that of radical spins by about hundred 
times hf at liquid Nitrogen temperature, the inter-system 
crossing is comparable with the relaxation of triplet, and 
the decay from the triplet to the ground singlet state is 
very slow which is about ms. In the following we will 
vary these parameters and take a close look at the effect 
of these relaxation processes. Another important param¬ 
eter is the exchange coupling between radical and triplet 
which can be both calculated by density functional the¬ 
ory [12 and estimated through the broadening of EPR 
spectra which is set to be —10 mK. In the following we 
will use mK as the energy unit. So in this section, we 
will first concentrate on SRTS and DRTS respectively, 
where we can vary the different relaxation parameters 
as our model could predict the time evolution of density 
operator and EPR spectra with any parameters in prin¬ 
ciple; meanwhile we present the 3-dimensional figure for 
the time-resolved EPR spectra. And then we will study 
the correlations between radical spins in DRTS caused 
by the random magnetic field driven by the relaxation of 
triplet. 

A. Time evolution of density operator 

Zero-field splitting D and E can be measured from the 
precursors; D is about 20mK and E 3mK. As we have 12 
parameters left for MonoTEMPO, in order to keep our 
model as simple as possible, for the moment we assume 
that g r = 9t — 2.0 because the observation of the inter¬ 
action between radical spin and triplet is somewhat inde¬ 
pendent of spin-orbital coupling which is strongly related 


to g-factor, and we also reck on that 7 * = 72,73 = 0 , 
7 ? = 721 7a = 0 as for high temperature this is almost 
true, yf = 7 |, and yf = 72 i.e., the inter-system cross¬ 
ing rate and the decay rate are spin-independent for the 
moment. Having some simplifications for these parame¬ 
ters, we are going to use some reason values captured by 
the previous literature for the relaxation processes and 
inter-system crossing, and our calculated exchange inter¬ 
action to compute the spectra to illustrate the physics 
inside the system, and then we can compare it with the 
experiments, so we can have some physical insights from 
these comparisons. From these simulations, we can see 
that the EPR spectra at high temperature can be pre¬ 
dicted by using our programs once if we have know the 
details of the relaxation and inter-system crossing. We 
set the parameters as: J ~ —lOmK, y\ = 0.067mK, 
y\ = 67.0mK, yf = 0.035mK, and V is about 0.67mK. 
The magnetic field is large compared with the above pa¬ 
rameters which is about 200mK. Now we can first com¬ 
pute time evolution of density operator withe different 
inter-system crossing rates. 

Inter-system crossing rate dependency. In this part, 
we simulate the case where the laser is turned on from 
time 0 , then turned off completely at t\ — 8.0 ns, and 
the density operator is evolved till the time t-z — 4000 ns. 
We first plot the triplet population as a function of inter- 
systern crossing rate yf as shown in Figure © because 
the population of triplet manifold is immediately related 
to the strength of EPR signals. We can see that the pop¬ 
ulation of triplet manifold will be higher when we tuned 
the inter-system crossing rate to be smaller. This is be¬ 
cause when the population transfers from singlet excited 
state to triplet, the zero-field splitting and the fast re¬ 
laxation of the triplet will transfer the population into 
Stotai = 2 state which will go back down to the ground 
much more slowly. So the more slowly the population is 
transferred, the more population will be left in the triplet 
manifold. 

Laser-held strength dependency. As shown in Fig¬ 
ure ©, we fix the inter-system crossing rate to be 33mK, 
and we see that as optical pumping is strengthened, there 
will be more population in the triplet manifold. 

Decay rate dependency. We can also show that the 
population of triplet manifold will drop faster when de¬ 
cay rate becomes larger as expected. The population of 
triplet manifold as a function of time with different decay 
rates is shown in Figure ©. 


B. EPR spectra 

As discussed in m we EPR spectra should be pro¬ 
portional to the response function of an open quantum 
system stated in equation 1211 From equation 1211 we can 
also see that EPR spectra is dependent both on time 
and magnetic field which results from the fact that the 
density operator is non-stationary for the unperturbed 
Hamiltonian. To simulate EPR spectra we take the per- 
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(b) 

FIG. 2: (Color online.) The population of the manifold with 
Soa = 1 as a function of time with varying inter-system cross¬ 
ing rate Rise = 7i for SRTS and Rise = 7i for DRTS is 
shown for both (a) SRTS and (b) DRTS. Notice that as Rise 
increases from 0.33 mK (black curve with square symbols), 
via 3.3 mK (red curve with triangle symbols) to 33 mK (green 
curve with circle symbols), the population of triplet becomes 
less during laser pumping for SRTS and DRTS. The slight 
difference between them is because there are more channels 
to dump the populations in triplet into singlet ground state. 


turbation B as the commutator [S x + s x , •] for SRTS and 
[5a, + Si 2 , + S 2 X , *] for DRTS, and the observable to be 
S x + $x for SRTS and S x + s± x + § 2 X for DRTS. Essen¬ 
tially, the physical meaning is that we put in a perturba¬ 
tion which flips the spins and then after some time we see 
the response of the system under the interaction with the 
environments. EPR spectra is Fourier transformation of 
the time-dependent response to the energy domain. 

Exchange-interaction-dependency. The main purpose 
that we use EPR for radical-triplet system is to look at 
how strong the exchange interaction between radical and 
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(b) 

FIG. 3: (Color online.) The population of the manifold with 
Soa = 1 as a function of time with different values of V is 
shown for both (a) SRTS and (b) DRTS; black curve with 
square symbols corresponds to V = 1.0 mK, red curve with 
triangle symbols V = 2.0 mK, and green curve with circle 
symbols V = 10.0 mK. Notice that as V increases, i.e., the 
laser pumping is strengthened, the population of triplet be¬ 
comes richer. 

triplet is. From Figure © we see that the broadening 
of EPR spectra is strongly dependent on the strength of 
exchange interaction J which is crucial both for quantum 
gate operations in QIP and molecular magnetic switch¬ 
ing in spintronics. In Figure © we keep all the other 
parameters the same as in Figure 0 with Rise = 33mK, 
V = 1.0, and to = 10 except varying J to OmK, 20mK, 
30mK, 50mK, lOOrnK, 200mK and 500mK. From this 
Figure, we can see that the line shape broadening of 
EPR spectra is strongly dependent on the strength of 
exchange interaction; as the exchange interaction be¬ 
comes stronger, the line shape will be broader. When 
J = 0 mK, the spectra is dominated by the perturbing 
oscillating field frequency. When J =£ 0, it will spread the 
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(a) 


(a) 




(b) 


(b) 


FIG. 4: (Color online.) The population of the spin mani¬ 
fold with Soa = 1 as a function of time with different values 
of decay rates Rd = 7 i for SRTS and Rd = 7 ® for DRTS 
is shown for both (a) SRTS and (b) DRTS; black curve with 
square symbols corresponds to = l.OmK, red curve with trian¬ 
gle symbols V = 2.0mK, and green curve with circle symbols 
V = lO.OmK. Notice that as Rd increases, the population of 
triplet decay faster. 


FIG. 5: (Color online.) The normalized EPR spectra as a 
function of magnetic field for both (a) SRTS and (b) DRTS 
are shown. The exchange interaction J is varied as taking 
values 0 (black square), 20 (red triangle), 30 (blue star), 50 
(green circle), 100 (magenta hollow circle), 200 (magenta hol¬ 
low triangle) and 500 mK (magenta hollow square). Notice 
spectra broadening and the width saturation. 


eigenvalues of the effective spin Hamiltonian, so we see 
the broadening. And when J ~ 500mK, the broadening 
will be saturated as the broadening is constrained by the 
spin-lattice relaxation time of triplet. 

Triplet relaxation-rate-dependency. In the previous 
calculations we always assume that the relaxation of 
triplet back to equilibrium is very fast. It is interest¬ 
ing to see what happens if we change the triplet relax¬ 
ation rate. Here we keep all the parameter the same as 
in Figure (0 except J = —50 mK and the triplet relax¬ 
ation rate Rt = taking the values 1.0 mK, 10.0 mK, 
50.0 mK, and 100 mK. As we see Figure 0, the relax¬ 
ation rate of triplet is also responsible for the broadening 
of EPR spectra. And the more important is that EPR 


spectra will be changed significantly if the triplet relax¬ 
ation is as slow as the relaxation of radical spin, i.e., 
yi ps yj. When Rt = 1 mK, the typical EPR spectra of 
a triplet with non-zero field splitting will appear. Mean¬ 
while, as the relaxation rate becomes larger, the triplet 
feature will be gone; we can only see the radical ^-spin. 

Analysis of EPR spectra by decomposing to the differ¬ 
ent spin correlations 

Time-dependent EPR spectra. The time-resolved EPR 
spectra is crucial for observations of the interactions 
between radical and triplet. Here we present the 3- 
dimensional picture for the calculated EPR spectra as 
a function of time and magnetic field. 




























V. 


CONCLUSION AND DISCUSSIONS 
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(a) 



(b) 

FIG. 6: (Color online.) The normalized EPR spectra of SRTS 
and DRTS as a function of magnetic field with varied triplet 
relaxation rate Rt taking values 1 (black square), 5 (red tri¬ 
angle), 10 (blue star), 50 mK (green circle). Notice spectra 
splitting due to the long-living triplet and its disappearance 
as the relaxation of triplet becomes faster. 
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( a ) 



(b) 

FIG. 7: (Color online.) The EPR spectra of SRTS as a func¬ 
tion of magnetic field with varied triplet relaxation rate R t 
taking values 67 mK (a) and (b) 0.67 mK while other param¬ 
eters are kept the same as in Figure 
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FIG. 8: (Color online.) The EPR spectra of SRTS and DRTS 
as a function of magnetic field with varied triplet relaxation 
rate Rt taking values 1 (black square), 5 (red triangle), 10 
(blue star), 50 mK (green circle). Notice spectra splitting 
due to the long-living triplet and its disappearance as the 
relaxation of triplet becomes faster. 
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